The recollapse problem of closed FRW models in 
higher-order gravity theories 



John Miritzis 

Department of Marine Sciences, University of the Aegean 
Electronic address: imyr@aegean.gr 

February 7, 2008 

Abstract 

We study the closed universe recollapse conjecture for positively curved FRW 
models with a perfect fluid matter source and a scalar field which arises in the 
conformal frame of the R + aR 2 theory. By including ordinary matter, we extend 
the analysis of a previous work. We analyze the structure of the resulted four- 
dimensional dynamical system with the methods of the center manifold theory 
and the normal form theory. It is shown that an initially expanding closed FRW 
universe, starting close to the Minkowski spacetime, cannot avoid recollapse. We 
discuss the posibility that potentials with a positive minimum may prevent the 
recollapse of closed universes. 

1 Introduction 



A closed Friedmann- Robert son- Walker (FRW) universe is often considered almost syn- 
onymous to a recollapsing universe. This is mainly due to our experience with the dust 
and radiation filled FRW models usually treated in textbooks. That this picture is mis- 
leading follows clearly from an example found by Barrow et al [I] according to which 
an expanding homogeneous and isotropic model with spatial topology S 3 satisfying the 
weak, the strong, the dominant energy conditions and the generic condition may not 
recollapse. Thus the problem of recollapse of a closed universe to a second singularity is 
delicate already in the FRW case. 
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The closed-universe recollapse conjecture states roughly that a closed universe cannot 
expand for ever, provided that the matter content satisfies some energy condition and 
has non-negative pressures. The conjecture was found true in certain spatially homoge- 
neous cosmologies j2], in certain spherically symmetric spacetimes j3] and in spacetimes 
admitting a constant mean-curvature foliation that possesses a maximal hypersurface [I]. 
In these investigations it has proved useful to demand that the dominant energy condi- 
tion and the positive pressure criterion hold (see also Ref. for a dynamical system 
approach) . 

In this paper we investigate the evolution of positively curved FRW models with a 
scalar field having the potential which arises in the conformal frame of the R + aR 2 
theory [HJ [7] and ordinary matter described by a perfect fluid with energy density p and 
pressure p. The motivation for this choice was presented in [S]. The purpose of the 
present article is to generalize the results in 8J by including ordinary matter and to 
correct the mistake found therein. 1 

The plan of the paper is as follows. In the next Section we write down the field 
equations, as a constrained five-dimensional dynamical system. We use the constraint 
equation to reduce the dimension of the system to four and after a suitable change 
of variables the system becomes quadratic. In Section III we analyze the structure of 
the equilibrium corresponding to the de Sitter solution using the methods of the center 
manifold theory. Furthermore, we find the so-called normal form of the dynamical system 
describing a large, slowly expanding universe with low total energy density; we show that 
such a universe cannot avoid recollapse. In the last Section, we consider potentials having 
a strict positive minimum and argue that this class of potentials prevent a closed universe 
from recollapse. 

2 Reduction to a 4-dimensional quadratic system 

In General Relativity the evolution of FRW models with a scalar field (ordinary matter 
is described by a perfect fluid with energy density p and pressure p) are governed by the 

1 In Ref. inequality (11) has the wrong direction (compare with 112[l in this paper). This mistake 
and a different rescalling (compare with JSJ) were the sources of the erroneous conclusion that an initially 
expanding universe avoids recollapse. In fact, inequality (11) must be reversed and as a consequence, 
the admissible trajectories of the system (16) start below the line H = \[2r in FIG.l. This implies that 
an initially expanded closed universe cannot avoid recollapse. Nevertheless, the calculations and the 
mathematical analysis of the system (12) near the equilibrium (0,0,0) remain correct. Moreover, the 
above mistake does not essentially affect the rest of the paper. 
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Friedmann equation, 

the Raychaudhuri equation, 

^ = -i(p + 3p + 20 2 -2v), (2) 
the equation of motion of the scalar field, 

+ 3-0 + V'(0) = 0, (3) 
a 

and the conservation equation, 

p + 3(p + p)- = 0. (4) 

a 

We adopt the metric and curvature conventions of jSj. a (t) is the scale factor, an 
overdot denotes differentiation with respect to time t, and units have been chosen so 
that c = 1 = 8nG. Here V ((f)) is the potential energy of the scalar field and V = dV/ dcj). 
We assume an equation of state of the form p = (7 — l)p, with 2/3 < 7 < 2. 
In what follows we assume that the potential function of the scalar field is 

V^) = V QO [l-e~^y (5) 

which arises in the conformal frame of the R + aR 2 theory [7j. 

The flat plateau of this potential is responsible for an early inflationary period of 
the universe and, for homogeneous spacetimes provides a mechanism of isotropization 
[3 EI]. From the field equations it is easy to see that in an expanding universe the energy 
density of the scalar field, namely 

E= l -^ + V(ct>) 

is a decreasing function of time. Since the energy density, p, of ordinary matter also 
decreases, it may happen that in a future time, E be comparable to p. In particular, for 
closed, k = 1, models, once the scale factor reaches its maximum value and recollapse 
commences i.e., H < 0, the term 3H<p in (jSJ) is no longer a damping factor, but acts as 
a driving force which forces the field to oscillate with larger and larger amplitude. If 
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Figure 1: The potential © 

this be the case, the repulsive effect of the cosmological term may drastically change the 
evolution of a classical FRW model. 

Setting <f) =: y, a/a =: H, we obtain from {I]) the constraint equation 

3H 2 + 3k/a 2 = p + ^y 2 + V(<f>), (6) 

which we use to eliminate a from the evolution equations Q-®- As a consequence, the 
dimension of the dynamical system is reduced to four and we obtain 

y = -3Hy-V'^), 
p = -3-fpH, 

H=\v^)~\y'-^^ P -H\ (7) 

We remind the reader that the exponential potential which is popular in the literature of 
scalar-field cosmologies has the nice property that V oc V, which allows the introduction 
of normalized variables according to the formalism of Wainwright et al [9 j. For an 
exponential potential the dimension of the dynamical system for a closed FRW model 
reduces to three [TT] . 

We simplify the system by rescaling the variables by the equations 

0^V3?20> y^V2V^y, p - ^p, H^^^H, t^Jj^t. (8) 
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Furthermore, in order to take account of the equilibrium point corresponding to the point 
at "infinity" and to remove the transcendental functions, it is convenient to introduce 
the variable u defined by 

u := e" , (9) 

and system (J7j) finally becomes 

u = -uy, 

y = — it + u 2 — 3Hy, 
p = -3>ypH, 

H = -(l-uf- -y 2 - ^^p - H 2 . (10) 

Note that under the transformation 0, the resulted four- dimensional dynamical system 
([TUjl is quadratic. 

Remark: The system (II U|) is not an arbitrary "free" four-dimensional system. In 
view of (jEJ) the initial conditions have to satisfy the condition 3Hq — p — \y 2 — V (<p ) < 0, 
or, in terms of the new variables, 

tfo-^-^/o-^(l-^o) 2 <0. (11) 

With a little manipulation of the equations (fTTH) it can be shown that, once we start with 
initial conditions satisfying at time t Q the inequality (|11|). the solutions of the system 
satisfy 

Hit) 2 -\ P {t)-\y(t) 2 -\{l-u{t)) 2 <0 

for all t > to. This is a general property of the Einstein equations, namely that the 
subsequent evolution of the system is such that the solutions respect the constraint. We 
conclude that the phase space of the system (fTUj) is the set 

S := | (u, y, p, H) G R 4 : H 2 - l -p - l -y 2 - 1 (1 - uf < J . (12) 

3 Stability analysis 

There are several equilibrium points of pOjl . Some of them correspond to static universes 
with a cosmological constant equal to \Aoo- in the study of the equilibrium points we 
note that u — 1 corresponds to <p — 0, i.e., to the minimum of the potential and u = 
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corresponds to = oo, i.e. to the flat plateau of the potential. In the following we pay 
attention to the most interesting equilibrium solutions which are: 

EQ1: (u — 0, y = 0, p = 0, H = 1/2) . This corresponds to the de Sitter universe with 
a cosmological constant equal to -\Aoo- We analyze the flow of (jl(Jj) near EQ1 in the next 
subsection. 

EQ2: (u — 1, y — 0, p — 0, H — 0) . This corresponds to the limiting state of an ever- 
expanding universe with H — > while the scalar field approaches the minimum of the 
potential and the scale factor goes to infinity. Equality in (|T2"|) which arises from the 
flat, k = 0, case defines a set on the boundary of E. We conclude that the point EQ2 
which corresponds to the Minkowski solution, is located on this boundary. The detailed 
structure of this equilibrium will be analyzed in subsection 13.21 

As we shall see, both equilibria are non-hyperbolic, i.e., some or all of the eigenvalues 
of the Jacobian have zero real parts. That means that the linearization theorem does not 
yield any information about the stability of the equilibria and therefore, more powerful 
methods are needed. The study of the qualitative behaviour of a dynamical system near 
a non-hyperbolic equilibrium point is difficult even in two dimensions. There are two 
general methods for simplifying a dynamical system having a non-hyperbolic equilibrium. 
The first is the center manifold theory. According to the center manifold theorem, 
the qualitative behaviour in a neighborhood of a nonhyperbolic equilibrium point q is 
determined by its behaviour on the center manifold near q. Since the dimension of 
the center manifold is generally smaller than the dimension of the dynamical system, 
this greatly simplifies the problem (cf. and also [13] for cosmological applications). 
The second method is the normal form theory, which consists in a nonlinear coordinate 
transformation that allows to simplify the nonlinear part of the system (cf. ^2] for a 
brief introduction). Both methods are used in the next two subsections. 

3.1 Center manifold for the system at EQ1 

It is easy to see that at the equilibrium point q = (u — 0, y — 0, p — 0, H — 1/2) , the 
Jacobian matrix of (|T0|) has one zero and three negative eigenvalues and, consequently 
the Hartman-Grobman theorem does not apply. The center manifold theorem implies 
that there exists a local 3-dimensional stable manifold through q (see for example |12j). 
That means that all trajectories asymptotically approaching q as t — > oo, lie on a 3- 
dimensional invariant manifold. Since q is a non-hyperbolic fixed point, the topology 
of the flow near q is non-trivial and is characterized by a one- dimensional local center 
manifold containing q. We prove the following result. 
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Proposition: The equilibrium point q = (0, 0, 0, 1/2) of (EJ) is locally asymptotically 
unstable. 

In order to determine the local center manifold of (fTUj) at q, we have to transform 
the system into a form suitable for the application of the center manifold theorem. The 
procedure is fairly systematic and will be accomplished in the following steps. 

1. The Jacobian of ffTOj) at q = (0, 0, 0, 1/2) has eigenvalues 0, —1, —3/2, and —37/2 

T T T T 

with corresponding eigenvectors (—2, 4/3, 0, 1) ,(0,0,0,1) ,(0,1,0,0) and (0,0, 3,1) . 
Let T be the matrix having as columns these eigenvectors. We shift the fixed point to 
(0, 0, 0, 0) by setting H = H — 1/2 and write (fTUj) in vector notation as 

z = Az + F(z), (13) 

where A is the linear part of the vector field and F (0) = 0. 

2. Using the matrix T which transforms the linear part of the vector field into Jordan 
canonical form, we define new variables, (x, 2/1, y 2 , 2/ 3 ) = x , by the equations 

u = — 2x, 
4 

y = 3 s + 1/2, 

P = 3y 3 , 

H = x + ?/i + ?/3, 

or in vector notation z = Tx, so that (J13)) becomes 

x = T^ATx + T -1 F (Tx) . 
Denoting the canonical form of A by B we finally obtain the system 

x = Bx + f(x), (14) 
where f (x) := T _1 F (Tx) . In components system (|14|) is 
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-|x z - xy 2 

fx 2 - y\ - \y\ + (3 7 - 1) y\ - 2x Vl - \xy 2 + (3 7 - 2) xy 3 + (3 7 - 2) 2/12/3 
fx 2 - \y\ - 4xyi - \xy 2 - Axy 3 - 3yty 2 - 32/22/3 
-37 (xy 3 + 2/12/3 + 2/1) 
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3. The system (JT5)l is written in diagonal form 

x = Cx + / (x,y) 

y = Py + g(x,y), (16) 

where (x, y) G R x R 3 , C is the zero lxl matrix, P is a 3 x 3 matrix with negative 
eigenvalues and /, g vanish at and have vanishing derivatives at 0. The center manifold 
theorem asserts that there exists a 1-dimensional invariant local center manifold W c (0) 
of (fTflJl tangent to the center subspace (the y = space) at 0. Moreover, W c (0) can be 
represented as 

r(0) = {(i,y)GlxM 3 :y = h (x) , \x\ < 6} ; h (0) = 0, Dh (0) = 0, 

for 5 sufficiently small (cf. ^2], P- 155). The restriction of (JTfij) to the center manifold 
is 

x = /(x,h(x)). (17) 

According to Theorem 3.2.2 in ^3], if the origin x = of ()17j) is stable (resp. unstable) 
then the origin of (JTfij) is also stable (resp. unstable). Therefore, we have to find the 
local center manifold, i.e., the problem reduces to the computation of h (x) . 

4. Substituting y = h (x) in the second component of (fTfij) and using the chain rule, 
y = Dh (x) x, one can show that the function h (x) that defines the local center manifold 
satisfies 

Dh (x) [f (x, h (x))} - Ph (x) - g (x, h (x)) = 0. (18) 

This condition allows for an approximation of h (x) by a Taylor series at x — 0. Since 
h (0) = and Dh (0) = 0, it is obvious that h (x) commences with quadratic terms. We 
substitute 





hi (x) 




d\X 2 


+ a 2 x 3 


+ 


(x 4 ) 


h{x) =: 


h 2 (x) 




bix 2 


+ b 2 x 3 


+ 


(x 4 ) 




. h (x) _ 




o 

C\X 


+ c 2 x 3 


+ 


(x 4 ) 



into (jTSjl and set the coefficients of like powers of x equal to zero to find the unknowns 

<2i,6i,Ci, .... 

5. Since y\ and are absent from the first of (|15|h we give only the result for h 2 (x) . 
We find b t = 32/27, b 2 = -32/81. Therefore, jnj) yields 

x = -^x 2 -|x 3 + |x 4 + 0(x 5 ). (19) 

It is obvious that the origin x = of (fT9*|) is asymptotically unstable (saddle point). The 
theorem mentioned after ()17j) implies that the origin x = of the full four-dimensional 
system is unstable. This completes the proof. 
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3.2 Normal form of the system near EQ2 

Regarding the stability of this equilibrium, it is easy to see that the eigenvalues of the 
Jacobian of (fTOj) are, ±z,0, 0, i.e., it is totally degenerate. Nevertheless, it is the most 
interesting case because in other equilibria the scalar field reaches the flat plateau, which 
is impossible if we restrict ourselves to initial values of H smaller than y/V^- We find the 
normal form of the system (jlU|) near the equilibrium point (u — 1, y — 0, p — 0, H — 0) . 
The idea of the normal form theory is the following: Given a dynamical system with 
equilibrium point at the origin, x = Ax + f (x) , where A is the Jordan form of the 
linear part and f (0) = 0, perform a non-linear transformation x — > x + h (x) , where 
h(x) = 0(|x| 2 ) as |x| — > 0, such that the system becomes "as simple as possible". 

To write the system in a form suitable for the application of the normal form theory, 
we shift the fixed point to (0, 0, 0, 0) by setting x = u — 1 and the system becomes 

± = -y - xy, 

y = x + x 2 — 3Hy, 

p = S-fpH, 

H=-x 2 - -y 2 - ^—^p - H 2 (20) 
4 T 6 H K J 

Now we make the non-linear change of variables 

2 3 7-2 3 

x -)• x - y T^P X + l H Vi 

lo 4 

3 7 - 2 3 TT 
y ^y + xy+ — py + -Hx, 

P, 

H^H + -xy, 

and keeping only terms up to second order, we obtain the normal form of the system, 
viz., 

3 

x = —y Hx, 

y 2 

y = x- -Hy, 

p = -3ypH, 

H = - 3 -^p- 1 -(x 2 + y 2 )-H 2 . (21) 
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Note that the results are valid only near the origin. 

Passing to cylindrical coordinates (x = r cos 9, y = r sin 9, p = p, H = H) , we have 

r = — rH, 
2 ' 

9 = 1, 

p = -37p#, 

H = -K 2 -H 2 -fc^p. (22) 
8 6 H K J 

We note that the 9 dependence of the vector field has been eliminated, so that we can 
study the system in the (r, p, H) space. The equation 9 = 1 means that the trajectory 
in the x — y plane spirals with angular velocity 1. The constraint (cf. (jl2}0 

h 2 <\p + \v 2 + \# 

becomes 

H 2 < i r 2 + -p. (23) 
We observe that the first and third of (|2*2*|) can be written as a differential equation 

d P o P 

ar r 

which has the general solution 

p = Cr 2 \ C > 0. (24) 
Therefore, for 7 = 1, we obtain from 



3 

r = — rH, 
2 ' 

H = --r 2 - -Cr 2 - H 2 , 
8 6 



It is convenient to rescale r by 



4C + 3 

so that the projection of (}2*2*j) on the r — H plane is 

r = — rH, 
2 ' 



24 

r, (25) 



H = -r 2 -H 2 . (26) 
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This system belongs to a family of systems studied in 1974 by Takens Note that 
the constraint becomes in the new variables 

H 2 < 2r 2 

and we conclude that the phase space of is given by 

-V2r < H < V2r. (27) 




H = -V2r 



Figure 2: The phase portrait of (J2EJ 

The phase portrait of (|2*Bj) is shown in FIG. 2 (see [H] for a detailed analysis). The 
system (I26|) has invariant lines H = cr with c = ±v2. Since no trajectory can cross 
the line if = cr, on any trajectory starting in the first quadrant below the line H = cr, 
H becomes zero at some time and the trajectory crosses vertically the r— axis. Once 
the trajectory enters the second quadrant, r increases and H decreases. At first sight, 
it seems probable that an initially expanding universe may avoid recollapse; in fact 
all trajectories starting above the line H = \/2r, asymptotically approaches the origin 
and the corresponding universes would be ever-expanding. But, (|2*Tj) implies that all 
trajectories with H > must start below the line H = y2r. In conclusion, inequality 
f!27|) leaves no room for an ever-expanding closed universe, contrary to what was claimed 
in 0. 

We conclude that for an initially expanding universe H continuously decreases while 
x and y oscillate with decreasing amplitude. H becomes zero at some time and the scale 
factor reaches a maximum value. Subsequently the universe begins to recollapse, i.e., H 
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continuously decreases below zero while x and y oscillate with increasing amplitude. A 
typical trajectory of (j21j) is shown in FIG. 3, where the variable p was suppressed. One 
obtains qualitatively the same picture for all 7 G [2/3, 2] . 



H 




Figure 3: A trajectory of (J2~T|) 



4 Further comments 

We have analyzed the qualitative behaviour of a positively curved FRW model filled 
with ordinary matter and containing a scalar field with the potential (J3J). This model is 
conformally equivalent to the positively curved FRW spacetime in the simplest higher 
order gravity theory, namely the R + aR 2 theory. We have shown that even for large 
initial values of H, near the flat plateau of the potential, the corresponding de Sitter 
equilibrium is asymptotically unstable. Furthermore, an initially expanding closed uni- 
verse in the neighbourhood of EQ2 cannot avoid recollapse. For open and flat models 
having potentials with a unique minimum, V (0) = 0, we have shown elsewhere ^H]; that 
in an expanding universe, the energy density p of ordinary matter, the Hubble function 
H and the scalar field (f) asymptotically approach zero. This theorem was proved without 
referring to the precise form of the potential. Putting all these results together, we may 
conjecture that potentials with a minimum equal to zero, cannot provide a mechanism of 
late accelerating expansion of the universe. On the other hand, in expanding universes 
with a potential having a positive minimum, the scalar field rolls down to the minimum 
of the potential and this residual cosmological term may explain the late accelerating 
expansion of the universe [T7j . 
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We illustrate this idea by the example of the more general quadratic theory, derived 
from the Lagrangian density R + aR 2 — 2A. The corresponding potential in the Einstein 
frame is 

V K {^) = V 00 {l-e'^^j +Ae~ 2 V^t ( 28 ) 

For every A > 0, the functions V\ (0) have the same qualitative behaviour as (jSJ) but, 
have a positive minimum, say V min at some m > 0. Both V m i n and (fi m increase with 
increasing A. 

We consider again expanding closed FRW models. It is easy to see that when V\ (0) = 
Vmin, the system (|7J) has an equilibrium ^0 = m , y = 0, p = 0, H = y^V mill /3j repre- 
senting the de Sitter solution. It can be shown simply by the linearization theorem 
that this equilibrium is stable. To avoid complicating expressions for the eigenvalues we 
proceed as in Section II and obtain the following system (compare to (fT0|) 'l 

u = -uy, 

y = —u (1 — u) + \u 2 — 3Hy, 
P = ~3lpH, 

H= l -(l-uf + -u 2 - -y 2 - - H 2 , (29) 

with A = A/Vqo. The equilibrium point 



1 + A' " ' r ' 2 V 1 + A 

corresponds to the de Sitter solution with a cosmological term equal to V m i n . Lineariza- 
tion of (|29|) near p is sufficient to show that this point is a sink (all eigenvalues of the 
Jacobian matrix of (|29p have negative real parts). Therefore p attracts all nearby so- 
lutions and initially expanding closed universes enter a phase of accelerating expansion. 
This attracting property of the de Sitter solution for expanding models is well known 
from the cosmic no-hair conjecture and is not restricted only to isotropic cosmology. We 
conclude that A = in (f2*Kj) is a bifurcation value for closed models that recollapse or 
not. 

However, de Sitter universe is not a global attractor for ([29)1 . Numerical experiments 
show that for highly curved models, or models filled with an excess of ordinary mat- 
ter, there are solutions of ()29j1 which recollapse. Conditions to prevent the premature 
recollapse of closed models were given in [TSj . 
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Our results are based on the analysis of the behaviour of the dynamical system 
(I10J) near the equilibrium solutions. A rigorous proof of the closed universe recollapse 
conjecture may come from the investigation of the global structure of the solutions of (J2J)- 
PJ) with k = +1. The study of the same question for Bianchi-IX models is an interesting 
challenge for mathematical relativity. 
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